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Abstract 

In the absence of a global description of the experimental data on elastic and soft diffractive 
scattering from the hrst principles of QCD, model-independent analyses may provide useful 
phenomenological insights for the development of the theory in the soft sector. With that in 
mind, we present an empirical study on the energy dependence of the ratio X between the 
elastic and total cross sections; a quantity related to the evolution of the hadronic central 
opacity. The dataset comprises all the experimental information available on proton-proton 
and antiproton-proton scattering in the c.m energy interval 5 GeV - 8 TeV. Generalizing 
previous works, we discuss four model-independent analytical parameterizations for X, con¬ 
sisting of sigmoid functions composed with elementary functions of the energy and three 
distinct asymptotic scenarios: either the standard black disk limit or scenarios above or 
below that limit. Our two main conclusions are the following: (1) although consistent with 
the experimental data, the black disk does not represent an unique solution; (2) the data 
reductions favor a semi-transparent scenario, with asymptotic average value for the ratio X 
= 0.30 ± 0.12. In this case, within the uncertainty, the asymptotic regime may already be 
reached around 1000 TeV. We present a comparative study of the two scenarios, including 
predictions for the inelastic channel (diffraction dissociation) and the ratio associated with 
the total cross-section and the elastic slope. Details on the selection of our empirical ansatz 
for X and physical aspects related to a change of curvature in this quantity at 80 - 100 GeV, 
indicating the beginning of a saturation effect, are also presented and discussed. 
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1. Introduction 


Run 1 of the CERN-LHC on pp collisions has revealed important new features of the 
strong interactions at the highest energy region reached by accelerator machines (7 TeV - 8 
TeV in the c.m. system) and excitement grows with the start of Run 2 at 13 TeV. However, 
despite all the theoretical successes of the Standard Model, the soft strong interactions (small 
momentum transfer) still constitute a great challenge for QCD. Once related to large distance 
phenomena, they are not formally accessed by perturbative techniques and the crucial point 
concerns the absence of a nonperturbative framework able to provide from the first principles 
of QCD a global description of the soft scattering states, namely the experimental data on 
elastic and inelastic diffractive processes (single and double dissociation) and that constitutes 
a long-standing problem [l|, Q . 

A renewed interest in these processes has been due to the large amount of experimen¬ 
tal data on elastic scattering (including total cross-section) that has been provided by the 
TOTEM Collaboration at 7 and 8 TeVpI-l^, more recently by the ATLAS Collaboration at 
7 TeV and also due to the present expectations with Run 2 at 13 TeV [^. Beyond the 
accelerator energy region, cosmic-ray experiments constitute the only tool for investigating 
particle properties and their interactions. Studies on extensive air showers (EAS) allow 
the determination of the proton-air production cross-section and from this quantity, it is 


possible to estimate the proton-proton cross-section Ifl, 111] at energies above 50 TeV. How¬ 


ever, in practice, the interpretation of the EAS development depends on extrapolations from 
theoretical formalisms that have been tested only in the accelerator energy region, resulting 
in rather large theoretical uncertainties, mainly in the estimation of the total cross-section. 

In the absence of pure QCD descriptions, the theoretical investigation of the soft strong 
scattering is still restricted to a phenomenological context and therefore characterized by 
some intrinsic limitations. In fact, although partial descriptions of the bulk of the experi¬ 
mental data can be obtained in the context of phenomenological models (for reviews see, for 


example, 12l-ll7]L the efficiency of any representative approach depends on a constant feed¬ 


back of new and adjustable parameters, which in turn are dictated by the new experimental 
data. As a practical consequence, equivalent data descriptions may be obtained in different 
phenomenological contexts, associated, in general, to different physical pictures for the same 
phenomenon (compare, for example, 181-129]). resulting, nearly always, in a renewed open 
problem. 

At this stage, another important strategy is the development of model independent de¬ 
scriptions of the experimental data, looking for quantitative empirical results that may work 
as an effective bridge for further developments of the QCD in the soft scattering sector and 
even selecting phenomenological pictures. As discussed in what follows, that is the point we 
are interested in here. 

Among the physical quantities characterizing the elastic hadron-hadron scattering, the 
ratio X between the elastic (integrated) and total cross sections as a function of the c.m. 
energy 


A'(s) = ^(s), 

^tot 


( 1 ) 
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plays an important role for several reasons. 


1. In the impact parameter representation (6-space), it is connected with the opacity (or 
blackness) of the colliding hadrons; as we shall recall, it is proportional to the central 
opacity (6 = 0) in the cases of the gray/black disk or Ganssian prohles. 


2. Due to model-dependence involved in direct measurements of the inelastic cross section, 
the s-channel unitarity constitutes the less biased way to obtain aineu namely from 
a tot — (^ei- Therefore, the same is true for the ratio ainpi / (Jtnt = I — X which, in turn. 


can be associated with the inelasticity of the collision 30 


3. Empirical information on X(s) and 1 — X(s) at the highest s and as s —)■ oo, provides 
crucial information on the asymptotic properties of the hadronic interactions, namely 
information on how atot and aei reach simultaneously their respective unitarity bounds 
(which is one of the main prospects in the LHC forward physics program [^). Empirical 
information on asymptopia is also important in the construction and/or selection of 
phenomenological models, mainly those based or inspired in nonperturbative QCD. 


4. Through an approximated relation (to be discussed later), the ratio X can be connected 
with the ratio between the total cross-section and the elastic slope parameter (H), 


F(s) 


1 O'tot 

Ibvr B 


( 2 ) 


This ratio, extrapolated to cosmic-ray energies, plays an important role in the study 
of extensive air-showers, as will be discussed along the paper. Moreover, it gives 
also information on the connection between atot and B at the highest and asymptotic 
energies. 


Based on these introductory remarks, our goal here is to present an empirical analysis 
on the ratio X (s) and discuss the implications of the results along the aforementioned lines, 
with main focus on asymptotic scenarios (s —)■ oo). In this respect, previous results obtained 
with smaller sized datasets and treating particular aspects, have already been reported in 


3ll-l34|: we shall review these results later. 

In the present work, our empirical study on the above quantities is once more updated, 
developed and extended in several aspects. The dataset on X comprises pp and pp scattering 
above 5 GeV, including all the TOTEM data at 7 and 8 TeV and also the recent ATLAS 
datum at 7 TeV. Generalizing previous works, we introduce four analytical parameterizations 
for V(s), consisting of two sigmoid functions composed with two elementary functions of 
the energy. Arguments and procedures concerning the selection of these suitable empirical 
parameterizations are presented and discussed in detail. As in 33|, [3^, we investigate all the 
three possible asymptotic scenarios: either the standard black disk limit or scenarios above 
or bellow that limit. Depending on the case and variant considered, the parametrization 
has only three to five free fit parameters. Based on the s-channel unitarity, the results for 
X(s) are extended to the inelastic channel, with discussions on the dissociative processes 
(single, double diffraction), including the Pumplin bound. We also treat the connection 
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between X(s) and the ratio b^(s), Eq. ([2]), together with discussions on the applicability of 
the results in studies of extensive air showers in cosmic-ray experiments. 

On the basis of the dataset, empirical parameterizations, £t procedures and results, 
our main four conclusions are the following: (1) although consistent with the experimental 
data, the black disk limit does not represent an unique or exclusive solution; (2) with the 
asymptotic X-value as a free £t parameter, all the data reductions favor a scenario below 
the black disk, with estimated average result 0.30 ± 0.12; (3) within the uncertainties our 
selected parameterization and £t result indicate that asymptopia may already be reached 
around 10^ TeV; (4) a change from positive to negative curvature in X(s) is predicted at 80 
- 100 GeV, which suggests a change in the dynamics of the strong interactions associated 
with the beginning of a saturation effect. 

The article is organized as follows. Our analysis on the ratio X is developed through 
Sections El tola where we present: the dataset investigated and the arguments why to treat 
also limits below and above the black disk (Section[2]); a detailed discussion on the selection of 
our empirical parametrization and variants considered (Section [3]); the fit procedures and the 
£t results (Section ia- In Section Owe discuss in a comparative way all the data reductions 
and present our conclusions on the £t results, which lead to the selection of two asymptotic 
scenarios: semi-transparent and black. In Section [6], we display the predictions for some 
quantities associated with the inelastic channel and the ratio between the total cross-section 
and the elastic slope. In Section [71 we discuss some possible physical aspects related to the 
analytical structure of the empirical parameterizations, presenting also comments on a semi¬ 
transparent asymptotic scenario. A summary and our final conclusions are the contents of 
SectionIHl In two appendixes, we present some basic formulas and results referred to along the 
text (Appendix A) and a short review of our previous works on the subject (Appendix B). 


2. Experimental Data and Asymptotic Scenarios 

2.1. Dataset 

Our dataset comprises all the experimental information presently available on the ratio 
X from pp and pp scattering above 5 GeV and up to 8 TeV (42 points: 29 from pp and 


35[ |. The statistical and systematic uncertainties 


13 from pp) and have been collected from 
have been added in quadrature. 

As we will show, the pp data at the LHG region play a crucial role in our analysis and 
conclusions. These 6 points are displayed in Table EJ four points at 7 TeV, one point at 8 
TeV obtained by the TOTEM Gollaboration and one point at 7 TeV recently obtained by 
the ATLAS Gollaboration. For further reference, a plot of the above-mentioned dataset is 
presented in Fig. [TJ As illustration, we have also included an estimation of the ratio X at 57 
TeV, evaluated from the experimental information on atot and aimi obtained by the Pierre 


Auger Gollaboration |3q. This value, that did not take part of our data reductions, reads 


A(V^ = 57TeV) = 0.317°;^^, 


where statistical, systematic and Glauber uncertainties [3^ have been added in quadrature. 
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Table 1: Experimental data on atot, o'ei and the corresponding ratio X = Uei/otot at the LHC energy region. 


(TeV) 

(^tot 

(mb) 

^ el 

(mb) 

^elj ^tot 

Collaboration 

7 

98.3 ± 2.8 

24.8 ± 1.2 

0.252 ± 0.014 

TOTEM [3j 

7 

98.6 ± 2.2 

25.4 ± 1.1 

0.258 ± 0.013 

TOTEM 

7 

98.0 ± 2.5 

25.1 ± 1.1 

0.256 ± 0.013 

TOTEM 

7 

99.1 ± 4.1 

25.4 ± 1.1 

0.256 ± 0.015 

TOTEM 

7 

95.4 ± 1.4 

24.0 ± 0.6 

0.252 dz 0.004 

ATLAS [7j 

8 

101.7 ± 2.9 

27.1 ± 1.4 

0.266 ± 0.016 

TOTEM [6] 
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Figure 1: Experimental data on the ratio X = aei/cFtot used in this analysis, from pp and pp elastic scattering 
in the energy region 5 GeV < ^/s < 8 TeV (accelerator data). At ^/s = 57 TeV it is shown, as illustration, the 
ratio obtained from the estimations of Utot and Uinei in cosmic-ray experiment (Pierre Auger Observatory). 


2.2. Asymptotic Scenarios 

We present here some theoretical, phenomenological and empirical arguments that led us 
to investigate all the three possible asymptotic scenarios. Let us denote A the asymptotic 
value of the ratio X: 


lim X(s) = A. (3) 

S^OO 

In what follows we discuss and display some numerical values for A, which will be considered 
as typical of each scenario to be investigated in Section 01 through data reductions. 
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2.2.1. The Black Disk 


As recalled in Appendix A the black disk limit reads 


A= 

2 


This limit represents a standard phenomenological expectation. It is typical of eikonalized 
formalisms, as the traditional models by Chou and Yang , Bourrely, Soffer and Wu jss]- 
40|, the hybrid approach by Block and Halzen 411 an d a number of models that have been 
continuously rehned and developed (for example, jl8l- l23L l26l-l28|). 


2.2.2. Above the Black Disk 

We have the arguments that follows for investigating scenarios above (beyond) the black 
disk. 


1. In the formal context, the s-channel unitarity. 


imposes 


^el I ^inel _ ^ 

^tot ^tot 

an obvious maximum bound for Y(s), namely 

A = 1. 


(4) 


2. In the phenomenological context, the [/-matrix unitarization scheme by Troshin and 


Tyurin predicts an asymptotic limit beyond the black disk ^ |43|]; in this approach. 


the reflective scattering mode becomes active at small values of the impact parameter 


44 


3. In the formal context, two well known bounds have been established for the total cross- 


section and the inelastic cross-section ^ 


TT 


TT 


^tot{s) < — In (s/so) and (Tinei{s) < (s/sq). 


mt 


Therefore, in case of simultaneous saturation of both bounds as s ^ oo, it is possible 
that Uinei/crtot —t 1/4, which from unitarity, Eq. (jl]), implies in 

4 

It should be noted that this limit does not correspond to an usual interpretation of the 
aforementioned asymptotic bounds. In fact, two fractions 1/2 are usually associated 
with each bound, in place of 1/4 in the inelastic case, which favors the black disk 
471. |49| . Even if a simultaneous saturation of both bounds might be 


scenario 


questionable j^, the number 0.75 can be considered as an instrumental choice for data 
reductions, lying between the black disk and the maximum value allowed by unitarity. 
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2.2.3. Below the Black Disk 

Nonetheless seeming a rather unorthodox possibility, in what follows we recall some 
results which suggest and indicate asymptotic limits below the black disk. 


1. In the publications of the TOTEM Collaboration, the authors quote the prediction 
for the total cross-section obtained by the COMPETE Collaboration (in 2002), with 


parametrization RRPL2, energy cutoff at 5 GeV and given by 51 




COMPETE 
tot 


(s) = 42.6 s 


- 0.46 


-33.4 s 


- 0.545 


+ 35.5 + 0.307 ln2(s/29.1). 


(5) 


where all coefficients are in mb and s is in GeV^. In the same publications, a £t to the 
Cel data above 10 GeV by the TOTEM Gollaboration is also presented j^, 

= 11.7- 1.59 ln(s) + 0.134 ln2(s). ( 6 ) 


Using the above parameterizations, the ratio X for s —)■ oo reads 

A = 0.436, 


which suggests a scenario below the black disk. We shall return to parameterizations 
Eqs. ([5]) and (E]) in the next sections. 


2. We have recently developed several analyses of the experimental data on atot, the p 


parameter and Cei 52l-l54l| . including the TOTEM Gollaboration results at 7 and 8 


TeV. For our purposes, we recall that the parametrization for the total cross section is 
expressed by 


atot{,s) = Regge terms + a + (3 ln'^(s/s/i) 


where a, (3 and 7 are free parameters. Fits to data on atot and p (using derivative 
dispersion relations) from pp and pp scattering above 5 GeV, have led to statistically 
consistent solutions either with 7 = 2 (hxed) or 7 > 2 (free fit parameter). In both 
cases, extension of the parametrization to a^i data (same 7 value) allowed to extract 

value A. In all cases we have obtained A < 

lowest central value around 1/3 (see a summary of the results in ji^, hgure 10, where 
the above point obtained with the TOTEM and GOMPETE results is also displayed). 
For future use, as a typical input in data reductions, we shall consider the lowest value 
obtained in these analyses, which, within the uncertainties, reads 


1 / 2 , with 


the ratio X (s) and its asymptotic 


A = 0.3. 


3. Model-Independent Parametrization 

As we shall demonstrate in the next sections, our model-independent parameterizations 
for X(s) constitute useful, practical and efficient tools in both the description of the ex¬ 
perimental data and the extensions to other physical quantities. Given this important fact, 
we present in this section a detailed discussion on the choices and steps that led us to the 
construction of the empirical ansatz for V(s). 








3.1. Empirical and Analytical Arguments 

3.1.1. General Aspects 

First, from Fig. [H as the energy increases above 5 GeV, the X data decreases up to the 
CERN-ISR region 20 - 60 GeV), where they remain approximately constant and then 
begin to increase smoothly. From a strictly empirical point of view, this rise in the linear-log 
plot scale may suggest a parabolic parametrization in terms of Ins, that is, an increasing 
slope (positive curvature). However, as discussed in Section 2.2, in the formal context, the 
Unitarity Principle demands the obvious bound 1 for X{s) as s —)■ cxd and hnite values are 
also dictated by the Froissart Martin bounds and all phenomenological models and empirical 
analysis, independently of the scenario (black disk, above or below). Therefore, except in 
case of existence of an unexpected singular behavior at some finite value of the energy, the 
above facts indicate a constant finite asymptotic limit for the ratio V(s), that is, a smooth 
saturation effect as s —)■ cxd. That in turn, demands a change in the sign of the curvature at 
some finite value of the energy, so that X (s) goes asymptotically to a constant limit with a 
decreasing slope (negative curvature). 

The above empirical and formal arguments concerning X(s) suggest an analytical param¬ 
eterization related to a sigmoid function (“S-shaped” curves), in order to impose the change 
of the curvature sign and an asymptotic limit (a constant). However, there is one more step, 
because the description of the data at low energies and the correct curvature dictated by 
the experimental data, depend also on a suitable choice for the functional argument of the 
sigmoid function. That led us to hrst express the ratio as follows. Denoting S a sigmoid 
function and / its argument depending on s, we express X as a composite function, 

X{s) = AS{f), f = f{s), (7) 


where, from Eq. ([3j), 


lim S{f{s)) = 1. 

s—>-oo 

Therefore and still on empirical grounds, the point would be to test different forms for 
S{f) and /(s) through fits to the experimental data, looking for statistically consistent 
descriptions with an economical number of free parameters. Obviously the problem does not 
have a unique solution and it is not possible to test all the analytical possibilities. 

However, the use of different classes of sigmoid functions combined with different classes 
of phenomenologically based elementary functions is a way to take into account the intrinsic 
uncertainty in the choice of the complete parametrization. With that in mind, we shall 
select two forms for S{f), combined (each one) with two forms for /(s), as explained in what 
follows. 


3.1.2. Sigmoid Functions S{f) 

Several classes of sigmoid functions have applications in different scientific contexts and 
that includes the logistic, hyperbolic tangent, error function, algebraic ratios and many 
others (see for example 1^; we shall return to its applications in Section 7.1). Here 
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we consider two classes of sigmoid functions. One of them, already used in all our previous 
analyses j3ll-l34|. is the Hyperbolic Tangent, denoted by 


= tanh/ = 


l-exp{-2/} 


l + exp{-2/} 

In addition, we now consider a Logistic function, with notation 

1 


S\f) = 


l + exp{-/}’ 


( 8 ) 


(9) 


which can be connected with by translation and scaling transformation. 


3.1.3. Elementary Functions f{s) 

We shall express f{s) in terms of elementary functions of the standard soft variable 
denoted 


In (s/so) = V, 


( 10 ) 


where sq is a fixed energy scale to be discussed later. Different functions and different 
conditions have already been investigated in our previous analyses 31 3^- Here, in order 
to generalize our previous parametrization and extend our analysis, we first express / as a 
sum of two terms: a linear function of the standard variable v and a function g{y) which, 
dictated by the data reductions, can account for possible deviations from linearity: 


/(s) f{v) = a + I3v + g{v), 


( 11 ) 


where a, ft and 7 are real free fit parameters. 

Different tests, using distinct datasets (only pp or including pp) and different energy cut¬ 
offs, under the above mentioned criteria, led us to choose two typical forms in soft scattering 
for g{v), namely either a Power-Law, 

gPL{v) = v\ ( 12 ) 

where 5 is an (additional) free £t parameter, or a Logarithmic-Law, 

gLL{v)=\\iv, (13) 


therefore, without additional parameter. 

Here ends the arguments and steps used in the analytical construction of our empiri¬ 
cal parameterizations (or ansatz) for W(s), which is summarized and denoted in the next 
subsection. 


3.2. Analytical Parameterizations and Notation 

Collecting and summarizing Eqs. ([7]) to flT^ and as a matter of notation, we shall 
express our four choices of parameterizations as follows. For each sigmoid function (logistic 
or hyperbolic tangent) we consider two variants associated and denoted by the functions 
9{v)- 
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- Logistic (L) with variant Power-Law {PL) 




A 


1 exp{-[a -h /31n(s/so) + 7 ln'^(s/so)]} ’ 
Logistic (L) with variant Logarithmic-Law (LL) 


(14) 




A 


1 + exp{-[Q; /91n(s/so) + 7lnln(s/so)]} 
Hyperbolic Tangent (PT) with variant Power-Law (PL) 


(15) 


-^PL (-s) = ^ tanh{a -1- /31n(s/so) + yhr (s/sq)}; 
- Hyperbolic Tangent (PT) with variant Logarithmic-Law (LL) 


( 16 ) 




= A tanh{Q; -t- /91n(s/so) + 7 lnln(s/so)}; 


with the condition 


( 17 ) 


s > So- 


At last, and still as a matter of a short notation for discnssion and data redactions, in 
what follows we shall refer to the sigmoid fnnctions as logistic or tank and to the variants 
as PL or LL. 


3.3. The Energy Scale 

In principle, the energy scale sq conld be considered as a free £t parameter. However, 
that introdnces additional nonlinearity in the data redactions and additional correlations 
among all the adjnstable parameters, which are not easy to control. In this respect, onr 
strategy was to £x the scale, bat simnltaneonsly, to investigate its effect in the £t resnlts. 
Specihcally, as snmmarized in Appendix B in onr previous analysis we have assumed either 


sn = 1 GeV^ (an usual choice in phenomenology) (MI,132| or sq = 25 GeV^ (the energy cutoff) 


33 I. l34| . Here, we consider the reasonable and efficient choice of sq as the energy threshold 


for the scattering states (above the resonance region), namely 


So = 4mp, 

where rop is the proton mass (see also [Hil, section 4.2, for further discussion on this en¬ 
ergy scale). We will show (Section 5.5.1) that the selected results do not depend on the 
aforementioned choices. 
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3.4- Constrained and Unconstrained Fits 

In applying parameterizations flT^ . flT^ . ffT6|l and flT7|) to data reductions there are two 
possibilities to treat the asymptotic parameter A. 


(1) to £x it to an assumed value (as one of those displayed in Section 2.2), imposing therefore 

an asymptotic scenario; 

(2) to treat A as a free parameter, leading to the selection of an asymptotic scenario. 


In the former case, which we shall refer as constrained fit {A fixed), there are only four 
free parameters in the case of the variant PL and only three with the LL. In the latter case, 
refereed to as unconstrained fit {A free) we have five free parameters with variants PL and 
four in the LL case. Certainly, that represents an economical number of free parameters: 
note that individual hts to atot and aeu as those in Eqs. ([5]) and ([6]), demand 10 or more 
parameters for the ratio X. 

As commented before, our previous analysis have been developed only with the sigmoid 
tanh and particular cases of the variant PL. For future discussion these results are summa¬ 
rized in Appendix B In what follows we present the fit procedures and all the £t results, 
which will be discussed, in detail, in Section O 


4. Fit Procedures and Results 

4.1. Fit Procedures 


The data reductions have been performed with the objects of the class TMinuit of ROOT 


Framework, with confidence level fixed at 68 % ^ . We have employed the default MINUIT 


error analysis with the criteria of full convergence The error matrix provides the vari¬ 


ances and covariances associated with each free parameter, which are used in the analytic 
evaluation of the uncertainty regions associated with htted and predicted quantities (through 


standard error propagation procedures [59 


As tests of goodness of fit we shall consider the reduced chi-squared, yfi jv., where v is the 
number of degrees of freedom and the corresponding integrated probability, P{x‘^,i') js^. 
However, the presence of both statistic and systematic uncertainties in the experimental data 
puts some limitation in a formal interpretation of these tests. For that reason, our goal is not 
to compare or select fit procedures or fit results but only to check the statistical consistency 
of each data reduction in a rather quantitative way. 

Since the parametrization is non-linear in at least three parameters, different initial values 
have been tested in order to check the stability of the result. In this respect we developed 
the procedure that follows. 

For each sigmoid (logistic or tanh) and each variant {PL or LL), Fqs. flTTD to (ITTD . we 
first develop the constrained fits {A fixed). In this case we consider the five numerical values 
displayed in Section 2.2 as representative of the three scenarios investigated, as summarized 
in the following scheme: 

- above the black disk —)■ either A = 1 (maximum unitarity) or A = 0.75 (a possible “formal” 
result); 
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- the black disk A = 0.5; 


- below the black disk —)■ either A = 0.436 (the result from the TOTEM and COMPETE 
parameterizations, Eqs. (E]) and ([ 6 ])) or A = 0.3 (lowest value we have obtained in j52l-l54|b 


For each hxed A, different initial values have been tested for the other parameters (4 
with PL and 3 with LL), until reaching stable convergence and consistent statistical results. 

In a second step, using the values of the parameters from the constrained fit result as 
initial values (feedback), we have developed the unconstrained fits, with each A now as free 
parameter too. 

As already stated, by hxing A we impose an asymptotic limit and by letting A as a free 
fit parameter, we select an asymptotic scenario. The following diagram summarizes the cases 
investigated: 

. , f logistic . f PL „ f constrained (A hxed) 

sigmoid < , , variant < ^ ^ =^ht< . ■ ^ / k r \ 

y tanh ( LL y unconstrained (A tree) 

In what follows we display the ht results with the logistic function followed by those with 
the tanh. All these results will be discussed in Section |3 


4 . 2 . Fit Results with the Logistic 
4-2.1. Variant PL 

In the case of the constrained ht, for each hxed A we have four ht parameters. The 
results and statistical information on each ht are displayed in Tableland the corresponding 
curves, compared with the experimental data, in Fig. Ha). For clarity, the plotted curves 
correspond only to the central values of the free parameters (i.e. without the uncertainty 
regions). 

With these results for the parameters as initial values, including each A value, the un¬ 
constrained hts have been developed. The results are displayed in Table E] and show that 
all data reductions indicate almost the same goodness of ht (x^/zz and P(x^,z/)) and the 
same asymptotic central value, namely A = 0.292. Although, in each case, the values of the 
parameters a, (3, 7 and 5 may diher, once plotted together all curves corresponding to the 
central values of the parameters overlap. That is shown in Fig. H^), where the correspond¬ 
ing uncertainty region has been evaluated through error propagation from the ht parameters 
(Table [3]) within one standard deviation. For comparison, we have also displayed the ratio 
X(s) obtained through the TOTEM and COMPETE parameterizations, Eqs. ([5]) and ([ 6 ]). 

4 . 2 . 2 . Variant LL 

The same procedure has been developed with the variant LL. In this case we have only 
three parameters {A hxed) or four (A free). The constrained ht results are displayed in 
Table 0] and Figure the unconstrained case in Table [5l Once more, all ht results 

converged to the same solution in statistical grounds and in the values of all the parameters. 
The corresponding curve including the uncertainty region is shown in Figure Hb)- 

4 . 3 . Fit Results with the Hyperbolic Tangent 

The same procedures discussed for the logistic function have been applied in the case of 
the tanh, as follows. 
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Table 2: Fit results with the logistic function, variant PL (Eq. (fTTll i and constrained case {A fixed), z/ = 38 
(Fig. ISIa)). 


A 

hxed 

a 

/3 

7 

(5 



0.3 

125.5(1.5) 

0.328(22) 

-123.7(1.5) 

1.56(12)x10-2 

0.811 

0.790 

0.436 

169.12(10) 

0.1828(72) 

-168.59(10) 

6.62(30)x10-3 

0.882 

0.677 

0.5 

211.37(96) 

0.1627(63) 

-211.156(95) 

4.74(21)x10-3 

0.899 

0.647 

0.75 

69.517(80) 

0.1315(51) 

-70.027(79) 

1.148(52) X10-2 

0.932 

0.589 

1.0 

82.898(75) 

0.1195(46) 

-83.825(73) 

8.78(40) X10-3 

0.944 

0.568 


Table 3: Fit results with the logistic function, variant PL (Eq. (imi l and unconstrained case {A free), u = 
(Fig. HKb)). 


A 

initial 

A 

free 

a 

/S 

7 

5 


P{X^^) 

0.3 

0.292(33) 

125.6(1.5) 

0.35(11) 

-123.6(1.5) 

1.66(51)xl0-2 

0.831 

0.756 

0.436 

0.292(33) 

169.9(1.3) 

0.35(11) 

-167.9(1.3) 

1.22(38) X10-2 

0.831 

0.756 

0.5 

0.292(33) 

212 . 2 ( 1 . 2 ) 

0.35(11) 

- 210 . 2 ( 1 . 2 ) 

9.8(3.0)xl0-3 

0.831 

0.757 

0.75 

0.292(32) 

70.5(1.9) 

0.36(12) 

-68.5(1.9) 

2.98(91)xl0-2 

0.832 

0.755 

1.0 

0.292(32) 

70.5(1.9) 

0.36(12) 

-68.5(1.9) 

2.98(91) X10-2 

0.832 

0.755 



Figure 2: Fit results with the logistic function, variant PL (Eq. (fTUl l: (a) constrained case {A fixed, Table 
[2]); (b) unconstrained case {A free, Table [3]) with the corresponding uncertainty region and the result from 
the TOTEM and COMPETE parameterizations, Eqs. ([S]) and ([5]). Legend for data given in Figured] 


4-3.1. Variant PL 

The fit results with the variant PL are displayed in Table [6] (constrained case), Table [7| 
(unconstrained case) and Fig. S] 
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Table 4: Fit results with the logistic function, variant LL (Eq. (fTKll l and constrained case (A fixed), ^ = 39 
(Fig. Ela)). 


A 

hxed 

a 

/3 

7 



0.3 

1.90(13) 

0.324(22) 

-1.95(14) 

0.789 

0.824 

0.436 

0.534(79) 

0.182(10) 

-1.122(78) 

0.858 

0.720 

0.5 

0.221(71) 

0.1620(91) 

-1.005(69) 

0.875 

0.692 

0.75 

-0.503(58) 

0.1302(71) 

-0.813(56) 

0.907 

0.638 

1.0 

-0.922(53) 

0.1186(65) 

-0.742(51) 

0.918 

0.617 


Table 5: Fit results with the logistic function, variant LL (Eq. (IlSp l and unconstrained case (A free), z/ = 38 
(Fig. Elb)). 


A 

initial 

A 

free 

a 

/3 7 

X^jv 


0.3 

0.293(26) 

2.05(59) 

0.346(87) -2.07(50) 

0.808 

0.794 

0.436 

0.293(26) 

2.05(59) 

0.346(88) -2.07(50) 

0.808 

0.794 

0.5 

0.293(26) 

2.05(59) 

0.346(88) -2.07(50) 

0.808 

0.794 

0.75 

0.293(26) 

2.05(59) 

0.346(88) -2.07(50) 

0.808 

0.794 

1.0 

0.293(26) 

2.05(59) 

0.346(88) -2.07(50) 

0.808 

0.794 



Figure 3: Fit results with the logistic function, variant LL (Eq. (ITHll 'l: (a) constrained case (A fixed. Table 
S]); (b) unconstrained case (A free, Tabled with the corresponding uncertainty region and the result from 
the TOTEM and COMPETE parameterizations, Eqs. ([S]) and ([5]). Legend for data given in Eigure[T] 


4^.3.2. Variant LL 

The fit results with the variant LL are displayed in Table [8] (constrained case), Table [9] 
(unconstrained case) and Fig. O 
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Table 6: Fit results with the tanh, variant PL (Eq. (fTHll l and constrained case {A fixed), v = 38 (Fig. SDa)). 


A 

hxed 

a 


7 

5 



0.3 

125.23(24) 

0.13243(93) 

-123.95(24) 

6.31(48)x10-3 

0.818 

0.780 

0.436 

169.21(13) 

0.0588(35) 

-168.50(13) 

2.23(15)x10-3 

0.844 

0.740 

0.5 

211.44(10) 

0.0477(28) 

-210.83(10) 

1.380(96) X10-3 

0.853 

0.724 

0.75 

68.06(14) 

0.0283(16) 

-67.67(14) 

2.56(18)x10-3 

0.870 

0.697 

1.0 

83.50(10) 

0.0204(16) 

-83.22(10) 

1.51(10)x10-3 

0.875 

0.687 


Table 7: Fit results with the tanh, variant PL (Eq. (IT5)) 1 and unconstrained case {A free), v = 37 (Fig. 
Hb)). 


A 

initial 

A 

free 

a 

(3 

7 

5 

X^/v 

P{x\^) 

0.3 

0.31(10) 

125.18(43) 

0.12(11) 

-123.99(43) 

5.7(4.9)x10-3 

0.838 

0.746 

0.436 

0.31(10) 

169.45(43) 

0.12(11) 

-168.25(42) 

4.2(3.8)x10-3 

0.838 

0.746 

0.5 

0.312(49) 

211.73(18) 

0.118(52) 

-210.54(17) 

3.3(1.4)x10-3 

0.837 

0.746 

0.75 

0.31(10) 

68.42(48) 

0.12(11) 

-67.23(48) 

1.04(0.93) xlO-2 

0.838 

0.745 

1.0 

0.31(10) 

83.92(46) 

0.12(11) 

-82.73(45) 

8.5(7.4)x10-3 

0.838 

0.746 



Figure 4: Fit results with the tanh, variant PL (Eq. (fTHI) !: (a) constrained case {A fixed. Table |6]); (b) 
unconstrained case {A free, Table [7]) with the corresponding uncertainty region and the result from the 
TOTEM and COMPETE parameterizations, Eqs. ([S]) and ([5]). Legend for data given in Figure [TJ 

5. Discussion and Conclusions on the Fit Results 

In this section we treat the ht results obtained with both the logistic (Section 4.2) and the 
tanh (Section 4.3). In the hrst three subsections, we summarize and discuss in a comparative 
way the data reductions with: constrained/unconstrained hts (Section 5.1), logistic/tanh 
(Section 5.2) and variants PL/LL (Section 5.3); after that we treat the inflection point 
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Table 8: Fit results with the tanh, variant LL (Eq. (I17|) l and constrained case [A fixed), i/ = 39 (Fig. EKa)). 


A 

hxed 

a 

/S 

7 



0.3 

1.289(54) 

0.1317(90) 

-0.786(58) 

0.796 

0.813 

0.436 

0.718(25) 

0.0587(35) 

-0.358(25) 

0.822 

0.777 

0.5 

0.605(20) 

0.0476(28) 

-0.291(20) 

0.831 

0.763 

0.75 

0.381(12) 

0.0283(16) 

-0.174(12) 

0.847 

0.737 

1.0 

0.2811(89) 

0.0204(11) 

-0.1259(87) 

0.853 

0.729 


Table 9: Fit results with the tanh, variant LL (Eq. (fT7)l ) and unconstrained case [A free), z/ = 38 (Fig. 
Uh)). 


A 

initial 

A 

free 

a 

/3 7 

X^jv 

P{x\y) 

0.3 

0.312(35) 

1.19(25) 

0.117(37) -0.70(21) 

0.815 

0.783 

0.436 

0.312(48) 

1.19(35) 

0.117(51) -0.70(29) 

0.815 

0.783 

0.5 

0.312(35) 

1.19(25) 

0.117(36) -0.70(21) 

0.815 

0.783 

0.75 

0.312(35) 

1.19(25) 

0.117(37) -0.70(21) 

0.815 

0.783 

1.0 

0.312(35) 

1.19(25) 

0.117(37) -0.70(21) 

0.815 

0.783 



Figure 5: Fit results with the tanh, variant LL (Eq. (IT7l) l: (a) constrained case {A hxed, Table E]); (b) 
unconstrained case {A free. Table E]) with the corresponding uncertainty region and the result from the 
TOTEM and COMPETE parameterizations, Eqs. ([S]) and dH]). Legend for data given in Figure [TJ 


(Section 5.4), present our conclusions on the fit results (Section 5.5), which will led us to the 
selection of two asymptotic scenarios: semi-transparent and black limits. 

5.1. Constrained and Unconstrained Fits 

In the case of constrained hts [A hxed), the values of the parameters, statistical infor¬ 
mation and curves are presented in Tabled Fig. [2](a) (logistic - PL), Table 01 Fig. [21(a) 
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(logistic - LL), Table El Fig. |ll(a) (tanh - PL) and Table [HI Fig. El(a) (tanh - LL). As 
already commented, for clarity, the curves in the hgures correspond to the central values of 
the parameters. However, error propagation from the £t parameters (one standard devia¬ 
tion) lead to typical uncertainty regions as those displayed in part (b) of each figure (to be 
discussed later). 

From the Tables, for i/ = 38 or 39, the values of the ^ lie in a typical interval 0.79 - 0.94 
and those of the P{x^, v) in the corresponding interval 0.82 - 0.57, indicating, therefore, sta¬ 
tistically consistent fit results in all cases investigated. We also note that the corresponding 
curves are consistent with all the experimental data analyzed (mainly within the uncertainty 
regions, not shown in the figures). 

In the case of unconstrained fits {A free), the value of the parameters, statistical infor¬ 
mation and curves are presented in Table El Fig. [2](b) (logistic - PL), Table El Fig. El(b) 
(logistic - LL), Table m Fig. |ll(b) (tanh - PL), Table |9l Fig. Etb) (tanh - LL). 

Here we note a remarkable fact: all fit results converged to practically the same solution, 
specially in what concerns the asymptotic limit. These values are summarized in Table flOl 
Small differences in the value of the other fit parameters are discussed in what follows. 

Table 10: Summary of the asymptotic limits obtained in the unconstrained fits, {A free) with the logistic, 
tanh, variants PL, LL and the corresponding i/ and 'X^ jv. It is also shown v and jv in the case of the 
black disk {A = 0.5). 


sigmoid: 

logistic 

tanh 

variant: 

PL 

LL 

PL 

LL 

A free 

0.292(33) 

0.293(26) 

0.31(10) 

0.312(48) 

u 

37 

38 

37 

38 

H 

to 

0.832 

0.808 

0.838 

0.815 

A hxed 

0.5 

0.5 

0.5 

0.5 

V 

38 

39 

38 

39 


0.899 

0.875 

0.853 

0.831 


5.2. Sigmoid Functions: Logistic and Hyperbolic Tangent 

First, let us compare the constrained results in part (a) of Figs. [2] and El (logistic) with 
those in part (a) of Figs. ID and [5] (tanh). We note that, above the experimental data, the 
rise of W(s) with the logistic is faster than with the tanh and the differences increase as A 
increases (from 0.3 up to 1.0). For example, in the case of A = 1 (fixed), at ^/s = 10^ GeV, 
X K. 0.55 with the logistic and X k, 0.45 with the tanh in both variants, PL and LL. 

On the other hand, in the unconstrained cases (part (b) of the same hgures), these 
differences are negligible, even in the asymptotic region. The small differences in the central 
values of A (namely ~ 0.29 for the logistic and ~ 0.31 for the tanh) are also negligible within 
the uncertainties in these parameters (Table (TU]). 

5.3. Variants PL and LL 

In the constrained fits, the values of the free parameters differ substantially with variants 
PL and LL, in both cases: logistic (compare Tables [2] and E]) and tanh (compare Tables El 
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andlH]). Despite these differences, once plotted together, all curves overlap, as can be seen in 
parts (a) of Fig. |2] {PL) and Fig. [3] {LL) for the logistic and parts (a) of Fig. 0] {PL) and 
Fig. |5] {LL) for the tanh. 

In the unconstrained case and variant PL (Table [3] for the logistic and Table [7] for the 
tanh) all results lead to practically the same asymptotic values, A = 0.292 (logistic) and A 
= 0.31 (tanh), with some differences in the central values of the other parameters. On the 
other hand, with variant LL the central values of the parameters are all the same up to three 
figures with the logistic (Table |5]) and also with the tanh (Table |9]). The same is true for the 
goodness of fit in each case. 

5 . 4 . Inflection Point 

By construction our parameterizations present a change of curvature (an inflection point), 
determined by the root of the second derivative. In what concerns the ratio between elastic 
and total cross section, this inflection suggests a change in the dynamics of the collision 
process and that will be discussed in Section 7.1. For future reference, we display in Table 
[TT]the position in the c.m. energy of the inflection point for all fits developed. We note that 
this position increases with the value of the asymptotic ratio A^ but it is restricted to an 
interval y/s ~ 80 - 100 GeV. We shall return to this important feature in Section 7.1. 


Table 11: Inflection point: roots of the second derivative of X{^/s) in GeV for each fit developed. 

logistic tanh 


A 

PL 

LL 

PL 

LL 

free 

81.6 

81.2 

80.1 

80.0 

0.3 

82.9 

82.4 

78.7 

78.5 

0.436 

92.2 

92.0 

86.7 

86.6 

0.5 

93.7 

93.5 

87.8 

87.8 

0.75 

96.5 

96.1 

89.7 

89.6 

1.0 

97.4 

97.1 

90.3 

90.2 


5.5. Conclusions on the Fit Results 

In what follows, based on the above discussions, we present our conclusions on the fit 
results, which will lead us to select one variant {LL) and two scenarios (semi-transparent 
and black limits) for further discussions and studies. 

5.5.1. Main Conclusions 

As regards constrained fits, given the relative large uncertainties in the experimental data 
and the small differences in the values of fv and P{')(^, v) for z/ = 38 - 39, we understand 
that all the fit results are statistically consistent with the dataset and equally probable 
on statistical grounds, even in the extrema cases oi A = 0.3 and A = 1. In other words, 
despite the large differences in the extrapolated results (at the highest and asymptotic energy 
regions), the constrained fits does not allow to select an asymptotic scenario. That leads to 
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an important conseqnence: althongh consistent with the experimental data, the black disk 
does not represent an nniqne or dehnite solntion. 

With regards to unconstrained fits, independently of the sigmoid fnnction (logistic or 
tanh) and variant {PL or LL), the data rednctions lead to nniqne solntions within the 
nncertainties, indicating a scenario below the black disk, with central valne of A aronnd 0.29 
(logistic) and aronnd 0.31 (tanh), as snmmarized in Tab. (TDl Given the convergent character 
of these solntions, we consider the unconstrained fits as the preferred results of this analysis, 
namely we nnderstand that the data reductions favor a semi-transparent (or gray) scenario. 
The small differences among onr fonr asymptotic resnlts may be associated with a kind of 
“nncertainty” in the choices of S{f) and f{s). In this sense, once constitnting independent 
resnlts (central valnes and nncertainties), we can infer a global {g) asymptotic valne, for 
example, throngh the arithmetic mean and addition of the nncertainties in qnadratnre (Table 

m- 


Ag = 0.30 ±0.12. 


(18) 


It is important to notice that, within the nncertainties, this resnlt is consistent with onr 
previous analysis using the tanh, ^ = 0.5 (fixed) and energy scale at 25 GeV^ (the energy 


cutoff), namely A = 0.36(8) in 33j and A = 0.332(49) in 3^. In other words, the convergent 


result does not depend on the energy scale. Moreover, within the uncertainties, the above 
value is in plenty agreement with the limits obtained through individual hts to atot and a^i 
data, using different variants and procedures (see Fig. 10 in Ref. 54l|). 

Goncerning the variants PL and LL, given the goodness of fits, the uniqueness character 
of the convergent solutions (represented by the same central values of the parameters up 
to three figures) and the smaller number of free parameters, we select the variant LL as 
our most representative result. From Tab. (TUI we can also infer another mean value now 
restricted (r) to variant LL and the two sigmoid functions: 


Ar = 0.303 ±0.055. 


(19) 


Within the uncertainties, this result is also in plenty agreement with the aforementioned 
analyses. In what follows we shall focus our discussion, predictions and extension to other 
quantities to this variant. 


5.5.2. Selected Results and Scenarios 

Although the unconstrained hts have led to an unique semi-transparent (or gray) sce¬ 
nario, we have shown that the constrained hts with the black disk present also consistent 
descriptions of the experimental data. Based on the ubiquity character of the black disk 
limit in the phenomenological context, we shall consider also this case as a selected scenario 
for comparative predictions and discussions. 

Summarizing, in what follows we concentrate on the results obtained with variant LL, 
sigmoid logistic and tanh and two scenarios, for short referred to as gray {A free) and black 
{A = 0.5). For comparison, in Fig. [6] we display all these results for W(s) together with 
the corresponding uncertainty regions and with the energy region extended up to 10^° GeV. 
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Numerical predictions at some energies of interest for X (s) (and other ratios to be discussed 
in what follows) are shown in the third columns of Table [T2] (logistic, A = 0.293 and A = 
0.5) and Table [13] (tanh, A = 0.312 and A = 0.5). 



Figure 6: Selected fit results for X{s) with variant LL (Eq. (ITKll i. A free (gray), A fixed (black) and sigmoid: 
(a) logistic; (b) tanh. 


Table 12: Predictions for different ratios at the LHC energy region and beyond, with the logistic, variant 
LL, A = 0.293 (gray) and A = 0.5 (black). 


A 

x/i (TeV) 

^elj ^tot 


^diff/ ^tot 

^totj B 


2.76 

0.2408(37) 

0.7592(37) 

0.2592(37) 

12.11(19) 


8 

0.2575(49) 

0.7425(49) 

0.2425(49) 

12.94(25) 

0.293 

13 

0.2637(65) 

0.7363(65) 

0.2363(65) 

13.26(33) 


57 

0.278(13) 

0.722(13) 

0.222(13) 

13.95(64) 


95 

0.281(15) 

0.719(15) 

0.219(15) 

14.11(74) 


2.76 

0.2366(32) 

0.7634(32) 

0.2634(32) 

11.89(16) 


8 

0.2625(44) 

0.7375(44) 

0.2375(44) 

13.20(22) 

0.5 

13 

0.2750(49) 

0.7250(49) 

0.2250(49) 

13.82(25) 


57 

0.3141(64) 

0.6859(64) 

0.1859(64) 

15.79(32) 


95 

0.3276(69) 

0.6724(69) 

0.1724(69) 

16.46(34) 


These predictions deserve some comments. From Fig. El in the case of A free and taking 
into account the uncertainty regions, the extrapolations are consistent with our asymptotic 
value 0.303 ± 0.055 above ^/s ~ 10^ TeV. In other words the results suggest that the 
asymptotic region might already be reached around 10^ TeV. That, however, is not the case 
with the black disk for which typical asymptopia are predicted far beyond 10^° TeV. In 
what concerns Run 2, let us focus on the predictions for X(s) at 13 TeV. The four results 
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Table 13: Predictions for different ratios at the LHC energy region and beyond, with the tank, variant LL, 
A = 0.312 (gray) and A = 0.5 (black). 


A 

\/i (TeV) 

^el / ^tot 


^diff / ^tot 

^tot/B 


2.76 

0.2404(36) 

0.7596(36) 

0.2596(36) 

12.08(18) 


8 

0.2577(46) 

0.7423(46) 

0.2423(46) 

12.95(23) 

0.312 

13 

0.2647(58) 

0.7353(58) 

0.2353(58) 

13.30(29) 


57 

0.282(11) 

0.718(11) 

0.218(11) 

14.16(57) 


95 

0.286(13) 

0.714(13) 

0.214(13) 

14.39(67) 


2.76 

0.2381(33) 

0.7619(33) 

0.2619(33) 

11.97(17) 


8 

0.2612(42) 

0.7388(42) 

0.2388(42) 

13.13(21) 

0.5 

13 

0.2719(46) 

0.7281(46) 

0.2281(46) 

13.66(23) 


57 

0.3038(57) 

0.6962(57) 

0.1962(57) 

15.27(28) 


95 

0.3145(60) 

0.6855(60) 

0.1855(60) 

15.81(30) 


from Tables [T^] and [12] are schematically summarized in Figure [T] Although the expected 
measurements at this energy might not strictly discriminate between the two scenarios it 
may possible to have an indicative in favor of one of them. Even if that were not the case, 
these new experimental information will certainly be crucial for a better determination of 
the curvature beyond the inflection point. 



Figure 7: Preditions for the ratio X at 13 TeV, with variant LL, gray (A free) and black {A = 0.5) scenarios 
and sigmoids logistic (circles) and tanh (squares). 
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6. Extension to Other Quantities 

With the selected results, let us present extensions and predictions to some other quan¬ 
tities of interest, that did not take part in the data reductions. 

6.1. Inelastic Channel: Ratios and Diffractive Dissociation 

The ratio between the inelastic and total cross-section is directly obtained via unitarity, 
namely 1 - X(s). The results (analogous to those for X(s) in Fig. [6]) are displayed in Fig. 
[HI together with the uncertainty regions and the experimental data available. As expected, 
all cases present consistent descriptions of the experimental data. Numerical predictions at 
the energies of interest for the ratio o'inei/o'tot are shown in the fourth columns of Table [TH] 
(logistic, A = 0.293 and A = 0.5) and Table US] (tanh, A = 0.312 and A = 0.5). 



Figure 8: Predictions, via unitarity, for the ratio between the inelastic and total cross sections, with variant 
LL (Eq. (lT5l) l. A free (gray), A fixed (black) and sigmoids: (a) logistic; (b) tanh. 

In what concerns the asymptotic limit, 

lim ^ = 1-A, 

s^oo atot 

contrasting with the black disk (1/2), our global {g) and restricted (r) estimations, Eqs. flT8|) 
and flT^ . predict 


1-Ag = 0.70 ±0.12 and 1 - A,. = 0.697 ± 0.055. 


Beyond elastic scattering, the soft inelastic diffractive processes (single and double dis¬ 
sociation) play a fundamental role in the investigation of the hadronic interactions. An 
important formal result on the diffraction dissociation cross-section concerns the Pumplin 
upper bound [Ty, l77j |: 


tot tot 


( 20 ) 
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where cr^*// stands for the soft inelastic diffractive cross section (the snm of the single and 
donble dissociation cross sections). In this context, the black disk limit (1/2) may be asso¬ 
ciated with a combination of the soft processes, namely elastic and diffractive, giving room, 
therefore to a semi-transparent scenario. 

In this respect, Lipari and Lnsignoli have recentely discnssed the experimental data 
presently available on these processes, calling the attention to the possibility that the Pnmplin 


bonnd may already be reached at the LHC energy region 78|. The argnment is based on 


a combination of the measnrements by the TOTEM and ALICE Collaborations at 7 TeV, 
which indicates 


^diff 
O'tot 


n 24+0-05 


Oel + O'diff 
O'tot 


— n 4QC^+0-05 

— u.^yo_Q Qg, 


^diff 

O'el 


= 0.952 


+0.20 

- 0.245 


snggesting therefore that the Pnmplin bonnd is close to satnration. 

In case of satnration, namely the eqnality in the above eqnation, it is possible to estimate 
the ratio Odiff/otot at the LHC energies and beyond. The nnmerical predictions for this ratio 
are shown in the hfth colnmns of Table [12] (logistic, A = 0.293 and A = 0.5) and Table [13] 
(tanh, A = 0.312 and A = 0.5). Obvionsly, the asymptotic valne of this ratio is zero in the 
case of a black disk scenario strictly associated with the elastic channel. 

Moreover, nsing the Pnmplin bonnd and Unitarity, we can also infer an npper bonnd for 
ratio Odiff/Oinei, namely 

g'j — ^ ~ 

Oinel 1 ^('^) 

The cnrves corresponding to this bonnd in all cases treated in this Section are shown in 
Fig. [9] together with experimental data. Once more, contrasting with the asymptotic nnll 
limit in a black disk scenario, onr global and rednced estimations {Ag and Ar), lead to the 
predictions: 



Rg = 0.29 ±0A2 and = 0.283 ± 0.055 for s ^ cx). (22) 


6.2. Ratio Y Associated with Total Cross-Section and Elastic Slope 

In cosmic-ray stndies, the determination of the pp total cross-section from the proton- 
air prodnction cross-section is based on the Glanber formalism jo il]. In this context, the 
nncleon-nncleon impact parameter amplitnde (prohle fnnction) constitntes an important 
ingredient for the connection between hadron-hadron and hadron-nnclens scattering. This 
fnnction is typically parametrized by 


aj{s,bj) = 


_ [1 + p(-s)] Ototjs) b 2 /j 2 B(s)] 
dvr B{s) 


(23) 


where p, atot and B demand inpnts from models to complete the connection. However, models 
have been tested only in the accelerator energy region and in general are characterized by 
different physical pictnres and different predictions at higher energies. As a conseqnence, the 
extrapolations resnlt in large theoretical nncertainties, as clearly illnstrated by Ulrich et al. 
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Figure 9: Upper bounds for the ratio between the dissociative (single plus double) and inelastic cross sections 
from the selected results for the gray and black scenarios (logistic and tanh). Experimental information from 
pp scattering at 0.2, 0.9 and 1.8 TeV (white marks) 
marks) 84-86|. 


79148,111 and pp scattering at 0.9, 2.76 and 7 TeV (black 


From the above equation, any extrapolation is strongly dependent on the ratio atot/B^ 
namely the unknown correlation between atot and B in terms of energy. 

One way to overcome this difficulty is to estimate the ratio 


Y{s) = 


through its approximate relation with the ratio X(s), treated in Appendix A 

Y (s) 167rX(s). 


(24) 


(25) 


That has been the motivation of the analysis presented in 31 


“an almost model- 

independent parametrization for the above ratio may reduce the uncertainty band in the 
extrapolations from accelerator to cosmic-ray-energy regions.” 

The behavior of Y (s) extracted in this way and in all cases treated in this Section are 
shown in Fig. [101 fh® corresponding numerical predictions at the energies of interest are 
displayed in the sixth columns of Table [I2] (logistic, A = 0.293 and A = 0.5) and Table IT3] 
(tanh, A = 0.312 and A = 0.5). 


7. Further Comments 

In this section we present some additional comments on two topics: possible physical 
interpretations associated with our empirical parameterizations (Section 7.1) and aspects 
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Figure 10: Predictions for the ratio between the total cross sections and the elastic slope, with variant LL, 
(Eq. ((HI)), A free (gray), A hxed (black) and sigmoid: (a) logistic; (b) tanh. 


related to a semi-transparent (or gray) asymptotic scenario (Section 7.2). 


7.1. On Possible Physical Interpretations 

As we have shown, the sigmoid and elementary functions selected for parameterizing 
X (s), represent suitable and efficient choices on both statistical and phenomenological grounds. 
Moreover, a striking feature is associated with the small number of free parameter involved: 
from 3 to 5, depending on the procedure. Despite the strictly empirical character of our 
analysis and strategies, let us attempt to connect the results provided by our sigmoid ansatz 
to some reasonable physical concepts. The goal is not a direct interpretation of the results, 
but to discuss some ideas that could contribute with further investigation on the subject. 
This discussion is not related to a specihc asymptotic scenario but focus only on our choices 
for the sigmoid and elementary functions. 

We discuss two aspects: (1) the possible connection of W(s) with the contribution of 
effective partonic interactions through the opacity concept; (2) the relation of this process 
with a saturation effect associated to a change of curvature in X (^. 

Beyond technology and exact sciences, sigmoid functions j^, 561 have applicability in a 
great variety of areas, as biological, humanities and social sciences, including neural networks, 
language change, diffusion of an innovation and many others 6nl-l62| . Application in analytic 


unitarization schemes (eikonal/U-matrix) has been considered by Cudell, Predazzi, Selyugin 
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65 


64| and a recent application related to polarized gluon density is discussed by Bourrely 


This class of functions is generally associated with the Pearl-Verhulst logistic processes, 
in which the growth of a population is bounded and proportional to its size, as well as, to the 
difference between the size and its bound. The process may be represented by the logistic 


differential equation 60|, 


dN 

dt 


= rN 


N 

M 


( 26 ) 
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where N = N{t) is the population at a time t, M is its maximum value (the carrying 
capacity) and r > 0 is the intrinsic rate of growth (the growth rate per capita). If r = r{t) 
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and/or M = M{t) the equation is with variable coefficients 

In our case, the logistic ansatz for S'(/), Eq. ([9]), is a trivial solution of the differential 
equation 


dS 

W 


= ^ [1 - . 


(27) 


In terms of the ratio X and variable v, Eqs. 
equation reads 

dv dv 


(0 and flTOl) respectively, this differential 


X' 

1-7 


(28) 


corresponding, therefore, to a logistic equation with variable coefficient df /dv. 

How to interpret this result in the context of elastic hadronic scattering at high energies? 
One conjecture is to look for possible connections with an effective number of partonic 
interactions taking part in the collision processes, as the energy increases. 

With that in mind, we recall that the wide class of QCD inspired models, presently in the 
literature, had its bases in previous concepts related to semi-hard QCD or mini jet models 
66l-l73|. The main ingredient concerned the separation of soft and semi-hard contributions 


in the scattering process; the former treated on phenomenological grounds and the latter 
determined under perturbative QCD, parton model and probabilistic arguments. 


Specifically, as recalled in Appendix A in terms of the opacity function (D), the proba¬ 
bility that an inelastic event takes place at b and s is given by 




-20 


(29) 


Denoting the probability that there are no soft (semi-hard) inelastic interaction by Pg {Psh), 
we can associate 

G.r^eiis, b) = l- PsPsH = 1 " (30) 

so that, the total opacity reads 


f!(s, b) = S)s(s, b) + {1 sh(s, b). 


(31) 


Let us now focus on the semi-hard opacity, which is constructed under probabilistic 
arguments and QCD parton model as follows. Let n{s, b) be the number of parton-parton 
collisions at s and b, which is associated with the probability of semi-hard inelastic scattering. 
Since Psh = is the probability that there are no semi-hard inelastic interactions, 

mean-free path arguments lead to the conclusion that the probability that hadrons do not 
undergo SH scattering can be expressed by 


P — 

nSH — e 


and therefore, from Eq. [301 the semi-hard opacity reads 

^SH{s,b) = ^n(s,6). 


(32) 
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In the theoretical context, n{s, b) is expressed in terms of parton-parton cross section 
and hadronic matter distribution (related to form factors in the g^-space). In general, the 
formalism involves a rather complex structure with different choices for the corresponding 
^s, ^SH, constituent cross-sections and form factors (see, for example, for a very recent 


version j. 

However, for our purposes, the main idea is that an effective number of partonic interac¬ 
tions is clearly connected with the opacity and the prohle functions and therefore, as recalled 
in Appendix A with the ratio X(s) (at least in what concerns the central opacity in gray 


disk and Gaussian prohles). In this context, the sigmoid behavior characterized by a change 
of curvature, suggests a change in the rate of effective partonic interactions taking part in 
the hadronic collisions. A fast rise at low energies is followed by a saturation effect starting 
at the inflexion point, which represents a change in the dynamics of the interaction. 

The energy ^/s of the inflection point in terms of the asymptotic ratio A and in all cases 
investigated is displayed in Tab. [HI We see that, despite the different values of A, the 
position of the inflection point lies in a rather restrict interval, namely y/s ~ 80 - 100 GeV. 

Goncerning this region, we recall that the UAl Gollaboration have reported the mea¬ 
surement of low transverse energy clusters (mini jets) in pp collisions at the GERN Gollider 
and y/s: 200 - 900 GeV [^. Extrapolation of the observed mini jet cross-section to lower 
energy (Fig. 10 in 75|) suggest that the region 80 - 100 GeV is consistent with the beginning 
of the mini jet production. Now, the rise of the mini jet cross-section has been associated 
with the observed faster rise of the inelastic and, consequently the total cross-section ( 75 


and references there in). Therefore, once X = 0^x1 a tot-, h seems reasonable to associate this 
behavior with a change of curvature in X and the beginning of a saturation effect. Our 
results are consistent with this conjecture. 

Although suggestive, we stress that the above arguments are certainly limited as effective 
physical interpretations of our analysis and results. However, we hope they may be useful 
as a possible context for further investigation. 


1.2. On a Semi-Transparent Asymptotic Scenario 

From the discussion in Section [5l our analysis favors the global result Ag = 0.3 ± 0.12; 
a scenario below the black disk. Let us discuss some results and aspects related to the 
possibility of a semi-transparent scenario. 

Up to our knowledge, in the mid-seventies, Fialkowski and Miettinen had already sug¬ 


gested a semi-transparent scenario in the context of a multi-channel approach ^ ; based on 


the Pumplin bound, this scenario has been also conjectured by Sukhatme and Henyey [88 


More recently, the advent of the LHG brought new expectations concerning the historical 
search for asymptopia, but interpretation of the data seems not easy in the phenomenological 
context. In respect to scenarios below the black disk in the pre and LHG era, we recall the 
facts that follows. 


1. Lipari and Lusignoli 89| and Achilli et ah 9^ have discussed the observed overestimation 


of a el (or underestimation of atnei) in the context of one channel eikonal models. That led 
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Grau et al. 2^ to re-interpret the Pumplin upper bound as an effective asymptotic limit, 


CTdiff 1 
^tot (^tot 2 


as 


s —)■ oo. 


Based on the behavior of the experimental data on X(s) and a combination of the results 
at 57 TeV by the Auger Collaboration together with the predictions by Block and Halzen, 
Grau et al. have conjectured a rational limit A = 1/3 as a possible asymptotic value 
therefore, consistent with our preferred results. 
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2. In the recent QCD motivated model by Kohara, Ferreira and Kodama 9ll-l94j the scatter¬ 
ing amplitude is constructed under both perturbative and non-perturbative QCD arguments 
(related to extensions of the Stochastic Vacuum Model). The model leads to consistent de¬ 
scriptions of the experimental data on pp and pp elastic scattering (forward quantities and 
differential cross sections) above 20 GeV, including the LHC energies and with extensions 
to the cosmic-ray energy region. The predictions for the asymptotic ratio X lie below 1/2 
and are close to 1/3; therefore in agreement with our selected results. 

3. As discussed before, combination of the parameterizations by the COMPETE {atot) 
and TOTEM {(Jei) Collaborations indicates the a^mptotic value A = 0.436. Moreover, 
combination of the ATLAS parametrization {(Jei) [Til with that by COMPETE reads A = 
0.456. Both, therefore, below the black disk. 

4. At last, we stress that through a completely different approach, several individual and si¬ 
multaneous fits to atot and p data, extended to fit the Uez data |52l-l54| , have led to asymptotic 
ratios in plenty agreement with our inferred global limit = 0.30 ± 0.12. 

All the above facts corroborate the results and conclusions presented here, indicating the 
semi-transparent limit as a possible asymptotic scenario for the hadronic interactions. 

At last, it may be instructive to recall a crucial aspect related to our subject. Besides 
empirical analysis on the ratio X, another important model-independent way to look for em¬ 
pirical information on the central opacity is through the inverse scattering problem, namely 
empirical fits to the differential cross section data and inversion of the Fourier transform 


connecting the scattering amplitude and the profile function ( Appendix A ). In this context, 
the Amaldi and Schubert analysis on the pp data from the CERN-ISR (23.5 < \fs < 62.5 


GeV) constitute a classical example of extraction of the inelastic overlap function 96 


Still restricted to pp scattering at the ISR, but including data from the Fermilab at 27.4 
GeV, several characteristics of the profile, inelastic and opacity function (in both impact 
parameter and momentum tran sfer space) have been extracted and discussed by Menon and 


collaborators in references 97H102| . The approach is characterized by analytical model- 


independent parameterizations for the scattering amplitude and analytical propagation of 
the uncertainties from the fit parameters to all the extracted quantities. An important 
observation from these analyses is the necessity of adequate experimental information in the 
region of large momentum transfer in order to obtain reliable results on the central region 
(small impact parameter). This crucial fact was already pointed out by R. Lombard in 
the first Blois Meeting on Elastic and Diffractive Scattering (1985): ”.. extrapolating the 
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measured differential cross section can be made in an infinite number of manners. Some 
extrapolated curves may look unphysical, but they cannot be excluded on mathematical 
grounds.’' 


95 


The experimental data from the ISR on pp collisions cover the region in momentum 
transfer up to 9.8 GeV^ and w ith t he inclusion of the data from the Fermilab, this region 

fallowing therefore consistent emp irica l analysis (see, for 


can be extended up to 14 GeV^ [lOO 
example. Fig. 7 in 98| or Fig. 1 in 


99[ and the detailed discussions in [100[). That, however. 


is not the case with the pp data (ISR, SppS Gollider or Tevatron) and, unfortunately, with 
the pp data from the LHG neither: the TOTEM measurements cover only the region up to 
2.5 or 3.0 GeV2. As a consequence, propagation of the uncertainties from the extrapolated 
hts led to rather large uncertainty regions at small values of the impact parameter. It is 
interesting to notice that this effect has been clearly pointed out in the theoretical papers 
by Khoze, Martin and Ryskin (see, for example. Fig. 2 in 


18| and related comments). 


It is also important to notice that even discrete Fourier transforms, through adequate bins 
intervals in momentum transfer, can not guarantee the consistency of the extracted prohle 
near the central region. Recall that one of the main controversy concerning differential cross 
sections at large momentum is related to the possibility of oscillation or smooth decrease 
and/or presence or not of a second dip. Obviously, all these important features (connected 
with the central region), are completely lost in empirical hts restricted to small or medium 
values of the momentum transfer (a fact that is implicit in the aforementioned statement by 
Lombard). 

On the bases of the above discussion, we un derst and that only experimental data in the 
deep-elastic scattering region (above ~ 4 GeV^ 1Q3|]) can provide the necessary information 


for reliable empirical extraction of the prohle function (and opacity) at and near the central 
region. In that sense, it would be very important if the experiments at the LHG could 
extend the regi on of momentum transfer. In this respect Kawasaki, Maehara and Yonezawa, 
stated in 2003 104j | (also quoted in 98|): “Such experiments will give much more valuable 
information for the dihraction interaction rather than to go to higher energies”. 


8. Summary, Conclusions and Final Remarks 

We have presented an empirical analysis on the energy dependence of the ratio X = 
aei/ctot, with focus on its asymptotic limit. The main ingredient concerned four analytical 
parameterizations constructed through composition of sigmoid and elementary functions of 
the energy. For each sigmoid, S (logistic or hyperbolic tangent), two elementary functions, 
/, were considered, given by a linear function plus a power law {PL) or a logarithmic law 
{LL) of the standard variable Ins/so, with sq = By expressing X(s) = AS{f{s)), 

with A the asymptotic limit, two ht procedures have been considered, either constrained 
{A hxed, imposing the asymptotic limit) or unconstrained {A as a free parameter, selecting 
the asymptotic limit). Based on empirical and formal arguments, we have investigated 5 
limits from A = 0.3 up to A = 1 (maximum unitarity), including the black disk case, A 
= 0.5. Altogether, we have developed 40 data reductions, associated with 5 asymptotic 
assumptions, 2 sigmoid functions composed with 2 elementary functions and 2 procedures 
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(unconstrained and constrained). The dataset comprised all the experimental data available 
on the ratio X above 5 GeV (and up to 8 TeV). This detailed empirical analysis led to two 
main results: 

Rl. Taking into account the statistical information and uncertainties, all constrained hts {A 
hxed) presented consistence with the experimental data analyzed, even in the extrema 
cases, A = 0.3 or A = 1.0. In other words, all the scenarios investigated are equally 
probable and our analysis does not allow to select a particular scenario (from the 
constrained hts). 

R2. Independently of the sigmoid or elementary function considered, all the unconstrained 
hts {A free) converged to an unique solution within the uncertainties, which is also 
consistent with the experimental data analyzed. Based on the four results, a global 
asymptotic value can be inferred: Ag = 0.30 ± 0.12; moreover, a restricted result to 
variant LL indicates A^ = 0.303 ± 0.055. 

These two results led us to two main conclusions: 

Cl. Although consistent with the experimental data analyzed, the black disk limit does 
not represent an unique or exclusive solution. 

C2. Our analysis favors a limit below the black disk, namely a gray or semi-transparent 
scenario, with A given above. 

Based on these two conclusions and given the ubiquity of the black disk scenario in eikonal 
models, we have presented along the text the results and several predictions considering both 
the black and semi-transparent cases with the logistic and hyperbolic tangent and restricted 
to our variant LL (since it presents the smaller number of free parameters). 

An important characteristic distinguishing the two scenarios concerns the approach to 
the asymptotic region. In the semi-transparent scenario asymptopia may already be reached 
around 10^ TeV and in the case of the black disk only far beyond 10^° TeV. 

In what concerns onr empirical parameterizations, it seems important to stress their 
efficiency in describing the experimental data analyzed, independently of the hxed physical 
valne of A and despite the noticeable small nnmber of free parameters. For example, the 
constrained black disk £t with the variant LL demands only three fit parameters and led 
to consistent description of all data. We stress, once more, the contrast with 10 or more 
parameters typical of individnal hts to atot and data. All that indicate the good quality 
of our analytical choices for S{f) and /(s), on empirical and phenomenological gronnds. 

Given the empirical efficiency of these analytical representations for X(s), we have at¬ 
tempted to look for possible physical connections with the underline theory/phenomenology 
of the soft strong interactions. Presently, as a first step, we can only devise some sugges¬ 
tive ideas relating a “population growth”, represented by the logistic differential equation, 
with the nnmber of effective parton interactions taking part in the collision. This number 
increases at low energies, tending to a satnration as the energy increases above the 80 - 100 
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GeV region. Despite suggestive, this conjecture is rather limited on physical grounds and 
must be investigated in more detail. 

It seems difficult to consider as accidental the semi-transparent scenario indicated in a 
vari^ of independent analyses as those presented here, as well as in our previous works 


33|, ISJ, l52l-l54| and other indications referred to in the text. 

It is worth mentioning that an asymptotic scenario distinct of the black disk has been 
recently d iscussed by I. Dremin, in the context of unitarity and a Gaussian prohle approach 


in5l - lin7l |. In this analysis the opacity is maximum at a hnite value of the impact parameter 


and the central opacity is determined by the parameter Z = AirB/atot- The asymptotic limit 
Z —)■ 0.5 corresponds to complete transparency, leading the author to propose a black torus 
(or a black ring) as a possible asymptotic scenario. Since, from Appendix A, Z{s) ~ l/4X(s), 
with our preferred restricted asymptotic result Ar (19), we predict 


0.83 


as 


s —)■ oo. 


corresponding, therefore, to a central semi-transparent core and not complete transparency. 
This asymptotic scenario, represented by a black ring surrounding a gray dis k, has bee n also 
predicted by Desgrolard, Jenkovszky and Struminsky, in a different context 108 . Il09| . 

At last, we stress that the new experimental data from Run 2 might not be decisive as a 
direct numerical selection of the scenarios discussed here (in terms of our predictions). How¬ 
ever, this experimental information on X{s) will be crucial for an improved determination 
of the curvature above the inflection point and consequently providing better accesses to the 
asymptotic region in model-independent analyses. 
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Appendix A. Basic Formulas and Results 

In this appendix we collect some formulas and results characteristic of the elastic hadron 
scattering, which are referred to along the text. We treat two subjects: the eikonal and 
impact parameter representations (Section A.l) and the approximate relation connecting 
the ratio between the elastic and total cross section with the ratio between the total cross 
section and the elastic slope (Section A.2). The text and notation are based on references 

QHEg. 

Appendix A.l. Impact Parameter and Eikonal Representations 
Appendix A. 1.1. Physical Quantities 

In elastic hadron scattering, the amplitude F is usually expressed as a function of the 
Mandesltam variables s and t = —q^. In terms of this amplitude the physical quantities of 
interest here, with the corresponding normalization, are the differential cross section 

0 = 7r|F(s,g)|2, (A.l) 
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the elastic (integrated) cross section, 


(^ei[s = 




the total cross section (optical theorem) 


(Jtot{s) = 47rJmF(s, g = 0), 


the inelastic cross section (via unitarity), 


(^ineli^^') ^ez('S) 


and the p parameter, 


p{s) = 


ReF{s, q = 0) 
Im F{s, q = 0) 


(A.2) 


(A.3) 


(A.4) 


(A.5) 


From the above formulas, the optical point reads 

d(^ ^ + P^] 

dq^ «2=n IStt 


(A.6) 


Appendix A. 1.2. Impact Parameter Representation 

The representation of the scattering amplitude in the impact parameter space is named 
profile function (denoted F) and in case of azimuthal symmetry, they are related by 


F{s,q) = i bdbJo{qb)T{s,b). 


(A.7) 


The unitarity principle in the impact parameter space is usually expressed in terms of 
the total, elastic and inelastic overlap functions, Gtot{s,b) = Gei{s,b) + Ginei{s,b), which, in 
terms of the prohle function reads 


2ReT{s,b) = \T{s,b)\^ + G,nei{s,b). 


(A.8) 


In this representation the total, elastic and inelastic cross sections are given, respectively by 


o'tot(s) = 4:71 / bdbReT{s,b), aei(s) = 27r / bdb\T{s,b)\‘^, 


(A.9) 


^inel ^TT / b db ^ . 


(A.10) 
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Appendix A.1.3. Eikonal Representation and the Opacity Function 
In the eikonal representation, the prohle fnnction is expressed by 

r(s,6) = (A.ll) 

where x(s, b) is the complex valned eikonal function. In this representation, from the nni- 
tarity relation Eq. flA.8j) . 

= (A.12) 

Since unitarity implies Imx{,s,b) > 0, we have Gin{s,b) <1, so that from Eq. flA.lOp . 
Gin{s,b) can be interpreted as the probability of an inelastic event to take place at given b 
and s. From this result, the imaginary part of the eikonal is associated with the absorption 
in the scattering process and for that reason (and the optical analogy), it is named opacity 
function, which we shall denote 


n(s,t) = Imx(s, 6). 


(A.13) 


Let us neglect the real part of the scattering amplitude, so that T{s,b) is a real valued 
function. In this case. 


r(s, 6) = 1 - (A.14) 

and by expanding the exponential term, in hrst order, 

T{s,b) ^n{s,b). (A.15) 


Therefore, the prohle is also connected with the hadronic opacity. In fact, in the optical 
analogy (Fraunhofer diffraction), the function p = = 1 — r(s,6) represents the 

modihcation of the incident wave caused by a diffracting object. In this context, without the 
object, p = 1 and T = 0 (no diffraction) and in case of a completely opaque object, p = 0 
and r = 1 (no transmission). 


Appendix A.l.f. Gray Disk, Black Disk and Gaussian Profiles 

From the above discussion, a real valued prohle function for a gray disk of radius R{s) 
and central opacity ro(s) = T{s,b = 0), is represented by 


From Eqs. 

by 


(jMl), 


^tot 


T{s,b) 


ro(s), 

0 , 


27rro(s)i?^(s) and aei 


0 <b < R{s) 
b > R{s) 

= 7iTI{s)R‘^{s), so that the ratio X is given 


Xgd — 


ro(g) 

2 


(A.16) 
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and the case of a black disk (ro(s) = 1) reads 


Xbd = 


(A.17) 


A Gaussian profile, with central opacity ro(s), 

is important because, through Eqs. flA.ip . flA.2p and flA.bp and flA.7p . it can describe the 
sharp forward peak in the differential cross section data, leading to the experimental deter¬ 
mination of the integrated elastic cross section j^, (see the next Subsection A.2). In this 
case, one obtains atot = 2vrro(s)i?^(s) and a^i = 7rrg(s)i?^(s)/2, so that 


^g(s) = 


ro(s) 


(A.18) 


Therefore, these examples provide an interpretation of the ratio X(s) as the energy 
dependence of the central opacity (or blackness) associated with the interacting hadrons. 


Appendix A.2. The X and Y Ratios 

The experimental data on the differential cross section is characterized by the dominance 
of a sharp forward peak, which can be well described by 


da 

dq^ 


da 

dq^ 



(A.19) 


with B the (constant) forward slope. Substituting in the optical point, Eq. flA.bjl . and then 
integrating Eq. flA.2j) . the elastic cross section reads 


(Tei{s) 


[1 + ^L(g) 

B{s) IbTT 


Since, from the experimental data p{s) < 0.14, to take 1 -|- 1 is a reasonable approxi¬ 

mation, leading to 


or, with our notation. 


B{s) ~ aUsV 


(A.20) 


Y (s) 167rX(s). 


(A.21) 


Appendix B. Short Review of Previous Results 

In this Appendix, using the notation dehned in Section [3l we review some previous results 
we have obtained in hts with the tanh, special cases of the variant PL and energy scales at 
1 and 25 GeV^. 
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In the 2012 analysis by Fagundes and Menon l32|, the dataset were restricted to pp 


scattering above 10 GeV and included only the hrst TOTEM datum at 7 TeV. In order to 
infer uncertainty regions in the extrapolation to higher energies, two extreme asymptotic 
limits have been tested by either hxing A = 1/2 (black-disk limit) oi A = 1 (maximum 
unitarity). The experimental data have been well described through variant PL with hxed 
5 = 2 , hxed Sq = 1 GeV^, and only three ht parameters: a, /3 and 7 (denoted 71 , 72 and 73 
Through the approximate relation Y ^ lOvrX (Appendix A.2), it was possible to 


m 


m) 


extend the extrapolation of the uncertainty regions to the ratio atot/Beu which, in the context 
of the Glauber model Ml], plays an important role in the determination of the proton-proton 
total cross-section from proton-air production cross-section in cosmic-ray experiments. 

In our subsequent analysis [33|, 13^ , the energy cutoff has been extended down to 5 GeV 


and the dataset included all the TOTEM measurements at 7 and 8 TeV 33(|, as well as 
the recent ATLAS datum at 7 TeV js^. Preliminary hts to only pp data with variant PL, 
energy scale as the energy cutoff (sq = 25 GeV^) and different A values led to an almost 
unique solution indicating the parameter 6 consistent with 0.5, within the uncertainties. We 
than hxed 6 = 0.5 and developed new hts including now the pp data, considering either 
constrained and unconstrained cases and diherent values for the A parameter. All data 
reductions presented consistent descriptions of the experimental data analyzed and in the 
case of the unconstrained ht we have obtained an unique solution with asymptotic limit 
below the black disk, namely A = 0.332 ± 0.049 34]■ 
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